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We present a direct comparison of the number of dynamically correlated molecules in the shear- 
mechanical and dielectric relaxations of the following seven supercooled organic liquids: triph- 
enylethylene, tetramethyl-tetraphenyl-trisiloxane, polyphenyl ether, perhydrosqualene, polybutadi- 
ene, decahydroisoquinoline, and tripropylene glycol. For each liquid we observe that the numbers 
of dynamically correlated molecules in the shear and in the dielectric relaxation are proportional. 
We show that this proportionality can be explained by the constancy of the decoupling index of the 
shear and dielectric relaxation times in conjunction with time-temperature superposition. Moreover 
the value of this proportionality constant is related to the difference in stretching of the shear and 
dielectric response functions. The most significant deviations from unity of this constant are found 
in a liquid with strong hydrogen bonds and in a polymer. 



INTRODUCTION 

The dynamical processes taking place in a supercooled 
liquid arc complex. This is due to the fact that the mo- 
tion of the molecules of a liquid close to the glass tran- 
sition is intrinsically collective. When the liquid enters 
this ultra- viscous regime pQ , the rearrangement of a par- 
ticle involves the motion of many of its neighbors. The 
idea that the dynamics becomes more and more cooper- 
ative has led to searches for a growing length scale as the 
dynamics slow-down upon cooling. So far no standard 
static correlation has revealed a detectable growing cor- 
relation length arising in the supercooled phase. On the 
other hand, dynamic correlations [2] may account for the 
evolution of the correlation length scales involved in the 
glass transition. Recently Berthier, Biroli and co-workers 
invented a simple and powerful method to estimate the 
four-point susceptibility, Xi [SHH]- The central idea is 
to estimate the four-point function via a more accessible 
three-point function. The function Xii which cannot be 
easily measured directly, can be approximated from the 
temperature evolution of any measured dynamic variable. 

We can measure the frequency-dependent dielectric 
susceptibility, e(w), and shear modulus, G(w) using the 
same cryostat and covering overlapping temperature- 
frequency ranges [SHE]- The piezo-shear-gauge (PSG) 
technique [S] allows us to measure the dynamic shear 
modulus of a supercooled liquid close to its glass transi- 
tion (where G typically assumes values between 0.1 MPa 
and 10 GPa) in a wide frequency range (10~ 3 — 10 4 Hz). 
In this work we extract and compare the number of dy- 
namically correlated molecules in the structural (alpha) 
relaxation from two different dynamic variables: the di- 
electric susceptibility and the shear modulus. 
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TABLE I: Properties of the seven liquids studied (from Refs. 
[91 4121 [TH] ). T g is the glass transition temperature, m is the 
Angell fragility, / is the temperature index (the values of I 
are reported for the highest and lowest temperature stud- 
ied) . The last column indicates if the liquid has or not a clear 
secondary /3-relaxation. All the data here refer to dielectric 
measurements. 



We analyze below dielectric and shear-dynamic 
data collected and published by our group [9l 
[TT) . available on-line [12] . This study focuses on 
seven liquids: triphenylethylene ("TPE"), tetramethyl- 
tctraphcnyl-trisiloxane ("DC704"), polyphenyl ether 
("PPE"), perhydrosqualene ("squalane"), polybutadi- 
ene ("PB20"), decahydroisoquinoline ("DHIQ"), and 
tripropylene glycol ("TPG"). DC704, TPE, PPE, 
squalane and DHIQ are molecular van der Waals bonded 
liquids, TPG has hydrogen bonds, and PB20 is a poly- 
mer with molecular weight of 5000 g/mol. All liquids 
were used as acquired. The PPE used is the Santovac®5 
vacuum pump fluid, and DC704 is the Dow Corning®704 
diffusion pump fluid. All the other liquids were acquired 
from Sigma-Aldrich. All the experimental details about 
these measurements can be found in Rcfs. [BHS]. 

Some properties of the liquids [5rfT2l ITS] are reported 
in Table [TJ Here the relaxation time r is defined by the 
inverse loss peak frequency and the glass transition tem- 
perature T g is defined as the temperature where the loss- 
peak is located at 27rl0~ 3 rad/s. The temperature de- 
pendence of the relaxation time around T g is expressed 
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in terms of the Angell fragility index [TBHTB] 



m = 



d\og 
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The temperature dependence of activation energy is 
quantified via the temperature index [T71 [TH] 



I(T) 



d1nAE(T) 
d\nT 



(2) 



where AE(T) is the activation energy defined by r = 
r exp(A£ , (T)/fcs7 1 ). Table [j] reports the variation of 
I(T) in the temperature interval studied [TS]. Table|l]also 
reports presence of a clear Johari-Goldstain /3-relaxation 
[19] in the dielectric spectrum of the liquid. 



COMPARISON OF SHEAR-MECHANICAL AND 
DIELECTRIC RESPONSES 

The studies carried out by our group [TTJ IS] 
focused on the temperature-dependence of the the shear- 
mechanical and dielectric a relaxation times (indicated 
with t q and r e , respectively). Furthermore, comparison 
of the shape of these two relaxation functions was pre- 
sented. The main conclusions of these studies may be 
summarized as follows: 

(i) The relaxation time of the shear modulus is gen- 
erally different from that of the dielectric susceptibility 
at the same temperature T. The shear-mechanical relax- 
ation is always slightly faster than the dielectric, r e (T) > 
tq(T). Nevertheless, the shear and dielectric charac- 
teristic alpha relaxation times evolve in a rather simi- 
lar way in the liquids studied when T is changed. This 
was discussed in detail in [5] where the decoupling index 
T e (T)/Tc(T) was reported and its insignificant tempera- 
ture dependence was established (t 6 (T)/tg(T) ~ const). 
This picture is also confirmed by other studies found in 
the literature I22H261. 



(it) The shear response function and the dielectric re- 
sponse function generally have different shapes. In liq- 
uids that do not show any detectable Johari-Goldstain 
/3-relaxation [19] the shape of each frequency-dependent 
response is found to be almost temperature independent. 
This feature is referred as time-temperature superposition 
(TTS), and it is found to hold to a very good degree in 
the temperature-frequency range explored 0120] both for 
the shear and the dielectric relaxation [31]. 

(Hi) For those liquids that have a clear beta-relaxation 
the alpha relaxations in the shear and dielectric spectrum 
seem to approach a temperature independent shape as 



the temperature is lowered. This has been presented in 
detail in Refs. [9l |20l [21j [31] suggesting that for the alpha 
process alone TTS applies, while in the full spectrum 
TTS is lost because of the presence of the beta process. 

In the following we show that (i) and (ii) imply that 
the shear and the dielectric numbers of dynamically cor- 
related molecules (for the liquids without clear beta re- 
laxation) are proportional in the temperature range stud- 
ied. Moreover, the same conclusion applies if we assume 
TTS (as suggested by (Hi) and also done in [35]) to hold 
for the alpha process in those liquids that have a sec- 
ondary relaxation. To understand this link it we first 
briefly recall how to approximate the four-point suscep- 
tibility. 



ESTIMATION OF THE NUMBER OF 
DYNAMICALLY CORRELATED MOLECULES 

The four-point correlator can be interpreted as the 
variance of the dynamics around its average value. One 
can estimate this function from the following equation 

(Refs. I3H3I2ZH2S]) 



Cp 



d\nT 



(3) 



In this equation yf is the normalized real part of the 
response function and cp is the configurational heat ca- 
pacity per molecule at constant pressure. The right-hand 
side of Eq. ^ is an approximation of Xi > it actually rep- 
resents a lower bound for this function. Nevertheless, 
this method is found to give values of the four-point sus- 
ceptibility in good agreement with the actual values of 
X4 when these can be evaluated directly (for example in 
computer simulations) [4j [5] . The characteristic value of 
the four-point function (i.e., the typical number of cor- 
related molecules in the relaxation, N cozr ) is associated 
with the maximum of %4, 



N COI r(T) = max[x 4 (w,T)] 



(4) 



The maximum of this function is consistently found to 
close to a frequency close to the loss-peak frequency of 
the alpha dynamics. 

The normalized response function x' appearing in Eq. 
([3]) is computed by subtracting a baseline parameter to 
the measured response (for example the dielectric suscep- 
tibility), subsequently dividing by the amplitude of the 
function l27l 



X(W,T) 



x(u,T) - Xc 
Ax 



(5) 
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FIG. 1: Real and imaginary part (top and bottom re- 
spectively) of the dielectric response of DC704 measured at 
T =229.6, 227.6 225.5, 223.5, 221.5, 219.5, 217.5, 215.4, 213.4 
K (open circles). The full lines are fits with the HN function 
Eq. |[6j|. The amplitude Ax and the baseline Xoo of the fitting 
function are shown graphically for the curve at T =229.6 K. 



correlated molecules can be obtained: the number of cor- 
related molecules in the shear relaxation Nq and in the 
dielectric relaxation N e . Note that for comparing these 
two numbers knowledge of cp is unnecessary, being only 
a constant multiplicative factor in Eq. Q. 

From Fig. [3]we can appreciate the growth of the shear 
and diellectric N co „ upon cooling. The minimum in- 
crease of N e is of a factor ~1.6 in TPE and its maximum 
increase is of a factor ~5.7 found in TPG. The relaxation 
times of the responses studied grow at least four orders 
of magnitude in all liquids. 

Let us now see the form assumed by the equations ^ 
and Q if TTS applies. To do this let us consider the very 
general expression for a (normalized) response function 
obeying TTS: 



x(w,T) = 0(«t(T)). (7) 

This is the case of Eq. [6] if the parameters a and /3 are 
kept constant. Differentiating the real part of Eq. ^ 
with respect to InT (as in Eq. (|3|) we obtain 



d(j)'{ujT(T)) dcj)'{x) dr{T) 



d\uT 



■■ uj- 



dx dlxiT' 



(8) 



where the prime indicates the real part and we have in- 
troduced x = ujt{T). The maximum of the function |8]) 
can be estimated setting w = t -1 (a minor correction 
term is present in the case of very large stretching [53]): 



In this work we fit all response functions with the 
Havriliak-Negami (HN) function [30] 



x(w,T) = Xoo + 



Ax 

[1 + (iTU>) a Y ' 



(6) 



The Appendix details how we introduce the assumption 
of TTS in the analysis of the spectra of the liquids pre- 
senting a beta process (this is done by fixing the a and 
/? parameters of the function ([6])). An example of the 
fitting is reported in Fig. [TJfor the dielectric responses of 
DC704. The normalized functions are shown in the upper 
parts of Fig. [2jA (dielectric) and in Fig. [2}B (shear). 
The function X4 obtained from Eq. [3] is shown in Fig. 
[2] for the dielectric responses and the shear-mechanical 
response of DC704. Note that the maximum of these 
functions at the same temperature is located at different 
frequencies in the shear and the dielectric case (as is also 
the case for the loss peaks of the responses). Moreover, 
the shape of X4 is slightly different in the shear and in 
the dielectric case as discussed in detail below. 

Once we have determined the maximum of X4> the 
quantity N COTI can be obtained via Eq. Q. In this way 
two independent estimates of the number of dynamically 



max 



aft'(W(r)) 
for 



/ 2 (i) 



d In r(T) 
c> InT 



(9) 

where f(x) = (d<j>' (x) / dx) . From Eq. (|9| it is clear that 
the growth of N is determined uniquely by the growth 
of the relaxation time upon cooling if TTS strictly holds. 
If (as stated in (i)) the decoupling index has a negligible 
temperature-dependence (r e (T) / 'tg(T) ~ const) then 



dlnT e {T) 
dlnT 



d In Tg (T) 
dlnT 



(10) 



This means that the decoupling index of the character- 
istic number of correlated molecules in the shear and di- 
electric relaxation is also constant as T is lowered: 



N e {T) 
N G (T) 



const 



(11) 



where the constant is determined by the stretching of 
the shear and dielectric relaxations. In other words, the 
growth of Nq and N e is identical upon cooling, while 
their difference in absolute values is set by the differ- 
ent (temperature independent) shape of the two response 
functions. 
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FIG. 2: (A- Top) Normalized real part of the dielectric response for DC704 at T = 229.6, 227.6 225.5, 223.5, 221.5, 219.5, 
217.5, 215.4 and 213.4 K. The full lines are the corresponding (normalized) fits to the HN functions. (A-Bottom) (cp/fcs)X4 = 
{dx'(u},T)/d\nT) 2 for the dielectric (left) and shear (right) relaxation at T= 228.6, 226.6, 224.5, 222.5, 220.5, 218.5, 216.4 
and 214.4 K. The full line are the corresponding (cp/fcs)x4 computed form the fitting HN functions. (B) Same as A for the 
shear-mechanical response 



COMPARISON OF THE SHEAR AND 
DIELECTRIC NUMBERS OF DYNAMICALLY 
CORRELATED MOLECULES 

The results expected from the constant decoupling in- 
dex and TTS can be readly checked. Plotting N e (T) 
versus Ng{T) as in Fig. [4] we check that these quantities 
approximatively differ only by a multiplicative factor in 
the temperature range studied for all liquids considered 
in this work. As seen in the inset of Fig. [4] all the data 
collapse onto the line N £ cx Nq if we multiply N e by 
the value {Nc/N e ) where the average is taken over the 
temperature range studied. 

We test Eq. (Ill directly in Fig. [5] As seen from this 
figure the decoupling index N e (T) /Nq (T) of the number 
of correlated molecules in the dielectric and shear relax- 
ation is very weakly temperature dependent and it does 
not show any clear trend of a systematic increase or de- 
crease. Note that, while N t (T)/Nc(T) stays constant, 
N £ and No both grow significantly upon cooling for all 
liquids (see Fig. [3}. 

In Fig. [6] we test further the constancy of the de- 
coupling index. Fig. [6] demonstrate the validity of Eq. 
(FLO}. The equality (91nr £ /91nT) 2 = (91nr G /91nT) 2 



seems to hold to a good approximation as expected from 
T e {T)/T G (T) ~ const. 

We stress once again that the multiplicative factor be- 
tween the shear and dielctric iV corr (i.e., the constant ap- 
pearing in Eq. (Ill) is determined by the shape of the 



shear and dielectric response. Indeed, the function /(l) 
defined in Eq. ^ depends on the form of the relaxation 
functions that set N e /N G ~ [/ e (l)// G (l)] 2 . When the 
response is modelled by the HN function ([6]), the value 
of /(l) depends only on the parameters a and j3: 



/(l) = -a0 , Re[i a /(l+i a ) 1+f> ] 



(12) 



From Eq. jl2| ) it is easy to understand that if the dielctric 
response function has approximatively the same shape 



as the shear-mechanical one, the constant of Eq. (Ill is 



close to unity. This is the case of DC704 as can be seen 
from Figs. [5]and[7] A for which Eq. ( 12 ) gives consistently 
N e /N G ~ [/ e (l)// G (l)] 2 ~ 1.35 using the values a and (5 
obtained from the fitting. 

If the dynamic shear modulus is instead much broader 
than the dielectric response (as shown in Fig. [7jB for 
TPG), N e /Nc is significantly larger than unity (as seen 
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FIG. 3: The quantity (cp/kB)N COII for shear (filled symbols) 
and dielectric (open symbols) as a function of (T — T g ) for 
the liquids studied (see legend). T g is the glass transition 
temperature for the dielectric relaxation from Ref. [5] (see 
legend and Table 

in Fig. [5]). Also for this liquid we can check that the 
value of this ratio is consistent with the equations given 
above obtaining N^/Nq ~ 1.97. 

Finally we want to stress that the stretching of the re- 
laxation function does not only significantly affect N COTT , 
but also the full shape of the function Xii 1 ^) as calcu- 
lated from Eq. This is illustrated in Fig. [8j If the 
stretching of the shear and dielectric response functions 
is similar as in DC704 (see Fig. [TjA), NJN G is close to 
one, but also the shapes of xli^) an d xl( w ) are quite 
similar as seen in Fig. [8jA. 

If the two response functions have significantly differ- 
ent stretching, the corrresponding four-point susceptibil- 
ities will have quite different shapes. This is the case of 
TPG whose dielctric response function is more "Debye- 
like" than the shear-mechanical one (Fig. [7]B). In this 
case the functions Xa is clearly broader than x% as seen 
in Fig. [8]B. 



FIG. 4: Characteristic number of correlated molecules in the 
dielectric relaxation plotted versus the characteristic number 
of correlated molecules in the shear mechanical relaxation at 
the same temperatures for the seven glass-forming liquids (see 
legend). The full line indicates N t = No while the dashed 
lines represents 7V e = XNg (with A =1.9 and 0.5 respectively 
for the upper and the lower line). In the inset we show the 
data collapse on the line N e oc Ng obtained when we multiply 
N € by the value (Ng/N € ) where the average is taken over the 
temperature range studied. 



less, these quantities are approximatively proportional in 
the explored temperature range. For five of the seven 
liquids studied the ratio between the shear and dielectric 
characteristic number of correlated moleculues is close 
to the unity. The most significant deviations from this 
unitary ratio are found in a liquid with strong hydrogen 
bonds and in a polymer. Finally, we showed that the dif- 
ference in these absolut numbers arises from the different 
stretching of the dielectric and shear response functions. 

Center for viscous liquid dynamics "Glass and Time" is 
sponsored by The Danish National Research Foundation 
(DNRF). 



CONCLUSIONS 

We have compared the shear-mechanical and dielec- 
tric characteristic number of dynamically correlated 
molecules for seven supercooled liquids close to the 
glass transition. The number of dynamically correlated 
molecules in the shear-mechanical relaxation is generally 
different from that of the dielectric relaxation. Neverthe- 
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laxation at the same temperature as a function of the (dielec- 
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APPENDIX: Analysis details 

Note that, although we consider the real part of x when 
computing the four-point susceptibility, we fit simultane- 
ously the real and imaginary part with Eq. ([6]). This is 
done minimizing the (generalized) residual \ 2 f° r a com- 
plex variable x, i.e.: X 2 = T,j( x Lp ~ x lh)* ' ( X L P ~ x lh) 
where the star indicates the complex conjugate. 
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Shear and dielectric response (Imaginary Part) of DC704 at T=219.5 K Shear and dielectric response (Imaginary Part) of TPG at T=198 K 
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FIG. 7: (A) Comparison of (normalized) imaginary part of the shear and dielectric response function of DC704 at the same 
temperature T =219.5 K (main panel). The amplitude of the functions is normalized to the loss maxima (e^ax an d Gmax) and 
they are plotted vs oj/ui ma x for an easier comparison of the shapes. In the inset we show the two susceptibilities as they are 
are measured. (B) Same as (A) for TPG at T =198 K. The black symbols indicate the last point included in the fitting (see 
Appendix). 



Shear and dielectric x 4 for DC704 at T=220.5 K Shear and dielctric x 4 for TPG at T=1 95 K 
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FIG. 8: (A) Comparison of the shear and dielectric four-point susceptibility of DC704 at the same temperature T =220.5 K 
(main panel). The amplitude of each function is normalized to the maximum (>T ax ) and it is plotted vs cj/cj max for an easier 
comparison of the shapes. In the inset we show x? an d X% calculated from Eq. Isl (B) same as (A) for TPG at T =195 K. 



A further remark on the computation of Xi is that Eq. 
(|3]) involves the derivative with respect to the tempera- 
ture that is, in practice, performed as finite difference. 
After obtaining the normalized curves we consider two 
successive functions measured at different temperatures 
T\ and T^{< T\). Each frequency scan of a response 
function is carried in a way that each curve has points 



in the same frequencies. The derivative appearing in Eq. 
([3]) is then computed, at the single frequency, as follows: 



dlnT 



In Ti — In To 



T=(T 1 +T 2 )/2 111 J l mJ 2 

Here we illustrate how we introduce the assumption 



of TTS in the analysis of the spectra of the liquids pre- 
senting a secondary relaxation process (as TPG, DHIQ, 
Squalane and PB20). This is done by fixing the stretch- 
ing parameters in ([6]) for the fitting in the following way. 

First we fit the spectrum at the lowest temperatures 
(with free stretching parameters) where the secondary 
process is well separated from the a relaxation. We ex- 
clude from this fit some of the high frequency data (af- 
fected by the secondary relaxation). To select which data 
to remove from the minimization we plot the logarithmic 
derivative of the imaginary part of the response (that is 
a = dlnx" I 'dhiuj) that shows a minimum a m i n at the 
frequency u> m i n where the a process meets the secondary 
relaxation. The frequencies larger than this uj m in are not 
considered in the fitting. 



Once the parameters are found from this low- 
temperature spectrum they are fixed to fit all the other 
spectra up to high temperature. In those fits we also ex- 
clude the high frequency points from the computation in 
the same way illustrated above. 

When this procedure is completed \i is estimated form 
the obtained fitting functions. All the liquids with a sec- 
ondary process have been treated in this way. Note that 
analyzing the data in this manner we are implicitly as- 
suming TTS. A final remark is that if we compute \i 
from the original (normalized) relaxation function, in- 
stead that from the fitting functions, we find a relative 
difference beteween the heights of the maxima only of a 
few percents. 



